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ABSTRACT 

Sufficiently massive clumps of molecular gas collapse under self-gravity and 
fragment to spawn a cluster of stars that have a range of masses. We investigate 
observationally the early stages of formation of a stellar cluster in a massive 
hlamentary infrared dark cloud, G28.34-1-0.06 PI, in the 1.3mm continuum and 
spectral line emission using the ALMA. Sensitive continuum data reveal further 
fragmentation in hve dusty cores at a resolution of several 10^ AU. Spectral line 
emission from 18o, CH3OH, ^^CS, H2CO and N 2 D+ are detected for the first time 
toward these dense cores. We found that three cores are chemically more evolved 
as compared with the other two; interestingly though, all of them are associated 
with collimated outflows as suggested by evidence from the CO, SiO, CH3OH, 
H2CO and SO emissions. The parsec-scale kinematics in NH3 exhibit velocity 
gradients along the hlament, consistent with accretion flows toward the clumps 
and cores. The moderate luminosity and the chemical signatures indicate that 
the five cores harbor low- to intermediate-mass protostars that likely become 
massive ones at the end of the accretion. Despite the fact that the mass limit 
reached by the la dust continuum sensitivity is 30 times lower than the thermal 
Jeans mass, there is a lack of a distributed low-mass protostellar population in 
the clump. Our observations indicate that in a protocluster, low-mass stars form 
at a later stage after the birth of more massive protostars. 
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Introduction 


Infrared dark clouds (IRDCs) are believed to represent the earliest stage in the for¬ 
mation of massive stars. Observations in different bands of the infrared wavelength have 
signihcantly improved our understanding of the dyna mics of these cRuds . The early identi- 
hcations of IRDCs based on the ISO and MSX data f Perault et ah 1996 : Egan et ah 1998 ) 
are extended by observations with the Spitzer Space Telescope. Higher resolution images in 
the Spitzer IRAC b ands reveal intricate de tails of IRDCs thanks to galactic plane surveys 
such as GLIMPSE flChurchwell et al.l 120091) . Observations by the Her schel Space Observa¬ 


tory provide a complementary view of I RDCs at far IR wavelengths (iMolinari et al.l 12011 


Wilcock et al.l 120121: iRaean et al.ll2012a|) . Some of the IRDCs remain in absorption at the 
70 /rm band of the Herschel Space Observatory, and transition to emission in the 250, 350 
and 500 /rm bands. These 70 /rm-dark IRDCs may represent dense clouds at the earliest 
evolution stage of star formation. 

Among the IRDCs, those with masses of 10^ Mq or higher at a scale of ~1 pc are the 
prime candidates to study massive star and cluster formation, which hereafter we refer to 
clumpj^. These clumps contain molecular masses similar to those t hat harbor massive 


as 


protostars JBeuther et a. 
HII regions lzhang et ah 


1998al: 


Kurtz 


2002 


Ketc 


> lO "^ Lq), the presence of complex molec ules flBeuther et a 


20111) . and massive and energetic outflows fjZhang et ah 


2002 a: iBeltran et ahll2004l: IWalsh et al.l 12014) and hypercompact 


3 2007), whose relatively large luminosities 


2001 


2007tJ: 


Rathborne et ah 2008 


Beuther et ah 


2002 bl) all point 


to active massive star formation. In contrast to clouds with embedded massive star forma¬ 
tion, IRDC clumps show few signs of star formation despite their large reservoir of molecular 
gas at high densities of > 10^ cm“^. A survey toward 144 IRDC clumps for H-;; 0 masers, a 


signp ost of star formation, found only 14% of maser occurrence in the sample fjWang et ah 


2006h. a frequency mu ch lower than that in high-mass protostellar objects and HII regions 
(iBeuther et ahl l2002al) . In addition, IRDCs have consistently lower gas temperatures and 
line widths. Studies in NH 3 found that they have a temperature on the o rder of 15 K, and 
a line width on the order of 2 km s“^ avera ged over a spatial scale of Ipc fjPillai et ahl 12006 


Wang et al.l 120081: iRagan et ahl I 2 OIII. l 2012 bl). lower than those in high-mass protostellar ob - 
jects (HMPOs) and HII regions I Molinari et ah 1998 : Beuther et ah 2002a : Lu et ah 2014). 
These characteristics, together with their low luminosities, place massive IRDC c lumps at 
an earlier evolutionary stage than HMPOs and HII regions flRathborne et ahl 120071) . 


^This paper follows the nomenclature used in IZhang et al.l (120091) . and refers a cloud as an entity of 
molecular gas of 10 — 100 pc, a molecular clump as an entity of < 1 pc that forms massive stars with a 
population of lower mass stars, and dense cores as an entity of 0.01 to 0.1 pc that forms one or a group of 
stars, and condensations as an entity of < 0.01 pc. 
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High resolution observations of IRDC clumps reveal structures at < O.lpc scales and 


tend to have high deuteration fractionation and CO dep 

etion ( 

Pillai et ah 2007. 

2012; 

Hernandez et ah 2011: 

Fontani et ah 

2011; 

Chen et ah 

2010|; 

Zhang et ah 

2 OO 9 I 1 . as expected 


20071 ). Thermodynamic analysis found that cores in some clumps are sub-virial, with 


their gas mass a factor of severa l greater than the virial mass de rived from the line width 


f Pillai et al.ll201ll: Ihi et ahl 120131: iTan et ahl 120131: ILu et ahl 12014 1. Considering that these 


cores are embedded in dense clumps with substantial external pressure, the imbalance be¬ 
tween the gravitational mass and the internal support can be even larger. These analyses do 
not include magnetic helds, which supply ad ditional support. Ind eed, in a recent polarization 


survey of 14 massive clumps with HMPOs, IZhang et ahl (1201411 found that magnetic fields 


fields derived from statistical analysis is t ypically 1 


2009 


10 mG ( 

Chart et ah 

2009; 

Tang et ah 

2014; 

Oiu et ah 

2014 

). The magnetic 


field strength required for a virial balance in these IRDCs cores, on the other hand, is on an 
order of 0.5 m G. Therefore, IRD C cores can still be in a virial equilibrium for a moderate 
magnetic field fjPillai et ahl 1201511 . 


Understanding the formation of massive cores is an essential precursor to improving our 
under standing of the fo rmation of star clusters. For a typical star-formation efficiency of 10 - 


30 % flLada et al.ll2010[l . and for a cluster having a Salpeter type initial stellar mass function 
(IMF), one would expect a 10^ Mq clump to form a cluster of mass between 100 to 300 Mq, 
with at least 75 stars of mass between 0.5 - 20 Mq. Therefore, significant fragmentation of 
molecular gas is required in the clump in order to form a cluster of stars. In such a cluster 
environment, the majority of stars are at stellar masses of < 1 Mq, which corresponds well to 


cm 


-3 


Bonnell & Bate 

2002 

Larson 

2005 

Teixeira et ah 

2007 

Zhang et al.ll2009 

Palau et ah 


20131) . The question arises when it comes to forming stars around 10 Mq or larger. Although 

oin T v"* /-I 1 v"* rc oo it r/~» oF t’o /^-PF /-w* T-* e~\ t r/~» v-v-i TO "X/T r^F e~t /-\'P ^ n 


based on direct high angular resolution observations fe.g. 

Cesaroni et ah 1999: Beltran et ah 

2006: Beuther et ah 2005: Zhang et ah 

2007; 

Jimenez-Serra et 

ah 

2 OI 2 II. thev mav not rep- 

resent the state of early fragmentation 

(Longmore et ah 

2011; 

Rathborne et ah 20081. Ob- 


servations of massive IRDCs provide direct insights into formation of massive cores. Indeed, 


significant fragmentation is observed in massive mo l ecular clumps (jZhang et al.ll2009l: ISwift 


20091 : IZhang fc Wang|l201ll: IWang et al.ll201ll. l2012l |2014J) . These studies find that massive 


cores are typically 10 to 10^ times more massive than the thermal Jeans mass, pointing to 
a signihcant support from turbulence as well as perhaps magnetic helds. In addition, cores 
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tend to further fragment, indicating that massive stars form in close groups. 


This paper extends the study of a massive IRDC G28.34-0.06 (hereafter G28.34) and 
presents sensitive ALMA observations of continuum and molecular line emission at the 230 
GHz band. The cloud, at a kinematic distance of ~4.8 kpc, contains sever al 10^ of dense 


gas along the infrared absorption hlarnent extending 6 p c in the sky flGarev et al.l 12000 


Rathborne et al.l 120061: iPillai et al.ll2006l: IWang et al.l 120081) . Figure 1 presents an overview 
of this region from the mid, far infrared to 1.3mm continuum, and gas in NH 3 . Two promi¬ 
nent dust continuum cl umps, PI and P2, are revealed in the 850 and 450 /xm images ob¬ 
tained from the JGMT flGarev et al.ll2000l) . Besides these two massive clumps, additional 
continuum peaks are revealed in the 1.2mm image obtained from the IRAM 30m telescope 


flRathborne et al.lI2006L120101) . Despite the fact that PI and P2 clumps contain a similar 


amount of dense gas within 0.3 pc, the P2 region has a high gas temperature of 30 K, large 
NH 3 FWHM line width of 3.3 km s“^ and is associated with a strong 24 /xm point source 



complex molecular chemistry. Furthermore, the gas in PI appears to be externally heated 
with temperatures decreasing from 20 K in the outside of the cloud to 13 K inside of the 
cloud. Likewise, the turbulence measured by the NH 3 line width s appears to decrease from 


pc scales to 0.1 scales. These observations led IWang et al.l (120081) to suggest that PI is at a 


much earlier stage of massive star formation compared to P2. 

The SMA observations at 1.3mm resolved the G28.34 PI clump into hve cores along 
the hlarnent with ma sses from 22 to 64 Mq and an average projected separation of 0.19 
pc fIZhang et al.ll2009l. hereafter Paper II). Subarcse cond-resolution ob servations in 0 . 88 mm 
reveal further fragmentation in two of the hve cores fjWang et al.l 1201 ll. hereafter Paper III). 
The masses of the fragments are typically a factor of 10 larger than the Jeans mass derived 
from the average temperature and density of the medium they are embedded in. Gollimated 
outhows are detected in the GO 3-2 emission, emanating from each of the hve co res whose 
temperatures indicate little heating by embedded protostars flWang et ahl l2012l. hereafter 
Paper IV). Besides GO, no m olecular line e missio n is detected in the 2-1-2 GHz band of 
the SMA at 230 and 345 GHz. Zhang et ah ( 2009 ) hnd strong GO depletion in the region. 
These previous studies found that cores in G28.34 PI clump likely harbors intermediate mass 
protostars, based on the energetics of molecular outhows. Furthermore, these cores follow 
turbulent fragmentation, rather than thermal fragmentation. 


Our previous observations of G28.34 with the SMA do not have adequate sensitivity 
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and/or dynamic range to reliably detect fragments of thermal Jeans mass (approximately 
2 Mq). In addition, the chemical properties in these cores are largely nnknown dne to a 
lack of detection of molecnlar line emission. Here, we present ALMA observations of G28.34 
PI in Band 6 (230 GHz). The continnnm data reach a la rms of 75 fiJy, or 0.065 Mq at 
1.3mm assnming a dust temperature of 15 K. There is a lack of a population of low-mass 
protostars with masses down to ~ 0.4 Mq within the clump. The observations also reveal a 
plethora of molecules arising from hve intermediate-mass cores. This paper focuses on the 
analysis of physical properties of the region, and potential population of low-mass protostars 
in the protocluster. The chemical properties of dense cores will be discussed in a future 
publication. The paper is organized as follows: Section 2 presents the observational setup. 
Section 3 presents results, and in Section 4, we discuss the implication of the data relevant 
to cluster formation. Section 5 presents the conclusion of the study. 


2. Observations 

G28.34 PI was observed with Atacama Large Millimeter/Submillimeter Array (ALMA) 
on 2012 November 15 in its extended conhguration with a total of 29 12-m antennas in 
the array (Project ID; 2011.0.00429.S). Including time for calibration, the observing session 
lasted about 1.5 hours. The actual time on-source is 34 mins. The precipitable water vapor 
(PWV) was about 1.5 mm, resulting in system temperatures of 80 to 100 K at the observing 
frequencies. The observations employed the Band 6 (230 GHz) receivers in dual polarization 
mode. The projected baselines range from 15m to 400m (11.5 - 308 kA). The resulting 
synthesized beam is 0.8" when using a robust weighting parameter of 0.5. The FWHM 
primary beam of the ALMA 12-m antennas is approximately 27" at the observing frequencies. 
To maximize continuum sensitivity, we used all 4 spectral windows with a bandwidth of 1.875 
GHz in each window. The correlator was set in the FMT mode which provides a uniform 
channel width of 488.3 kHz (0.62-0.67 km s“^ at the observing frequencies). The frequency 
coverage of the data range from 218.07 to 219.95 GHz in SPWO, from 215.62 to 217.49 
GHz in SPWl, from 230.18 to 232.06 GHz in SPW2, and from 232.82 to 234.70 GHz in 
SPW3, respectively. Quasars J1751-I-096 and J1830-I-063 were used for bandpass and time 
dependent gain calibration, respectively. The flux calibration was achieved by observations 
of Neptune. The visibility data were calibrated in GASA by the ALMA supporting staff. 
We construct continuum visibility data using the line free spectral channels for the science 
target G28.34 PI. The continuum emission was self calibrated in GASA to further improve 
the dynamical range in the maps. The gain solutions from the self calibration were applied to 
the spectral line data. We found that the improvement in dynamical range in the continuum 
is about 20%. The continuum image reaches a la rms noise of 0.075 mJy in a synthesized 
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beam of 0".85 x 0".64. The spectral line sensitivity per 0.7 km s“^ channel is 2.5 mJy for 
SPWO, SPWl and SPW3. The noise in SPW2 is a factor of 3 higher for unknown reasons. 
Thus we exclude the SPW2 data for the continuum image. 


Results 


3.1. Continuum Emission 


Figure 2 presents the 1.3mm continuum emission of G28.34 PI. The ALMA image 
reveals 5 continuum emission peaks along the infrared dark hlament. The peaks of the 
continuum emission coincide w ith the 1.3mm SMA continuum sources SMAl through 5 
reported bv IZhang et al.l (120091) . There appears to be another continuum peak 5" southwest 
of SMA5. This dust peak is located at 10% of the primary beam response of the ALMA 
12-m antenna. With a primary beam corrected p eak flux of 4.8 mJ y, it is below the 3cr 
sensitivity limit of the previous SMA observations (jZhang et al.l 120091) . Since this feature is 
at 10% of the primary beam response, clean does not work well. We refrain from discussion 
of this feature. 

The superior sensitivity of the ALMA data reveals additional features in the contin¬ 
uum emission. SMAl, a single cont inuum source in the 1.3mm and 0.8 mm SMA data 


fjZhang et al.l 120091: IWang et al.l 1201 ih . becomes a group of 4 continuum peaks. The bright¬ 
est one coincides with the SMAl position. The other 3 continuum peaks are at fluxes of 
4.0, 2.0, 1.9 mJy, respectively. In addition, spatially extended emission is detected around 
SMA2, SMA3 and SMA4. Although SMA5 is a single continuum peak in the SMA data, 
the ALMA image reveals another continuum peak 2" to the east of SMA5. Overall, ALMA 
image reveals fainter and spatially extended emission. The integrated fluxes in the ALMA 
continuum image is 228 mJy, which is about 60% h i gher t han the integrated flux of 143 mJy 
in the SMA 1.3mm image reported in IZhang et al.l (1200911 . 


For convenience of cross referencing, as well as following the convention in the previous 
papers, we name the 5 groups of sources in the AL MA continuu r n ima ge Core 1 through 
Core 5, corresponding to SMAl through SMA5 in IZhang et ahl (120091) . respectively. To 
identify structures in the 1.3mm continuum emission, we employed dendrogram algorithm 
to decompose the continuum emission. However, the decomposition misses some obvious 
emission features. We then resort to fitting elliptical Caussian to the continuum image, but 
using the smallest features identified by dendrogram as a starting point of initial guesses for 
the Caussian fitting. We identify a total of 38 components in the simultaneous Caussian fit. 
The rms in the residual image is 6 x 75/iJy, corresponding to a detection limit of about 0.4 
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Mq. We report the coordinates, flux, the size of the objects in Table 1. The association with 
previous reported SMA continuum peaks is also listed in the table for cross referencing. The 
brightest components with peak fluxes > 1.9 mJy beam“^ are also labeled in Figure 2. 


Wang et al.l (120121 1 derived a gas temperature of 9 to 22K using the higher angular 


resolution NH3 (1,1) and (2,2) data. Adopting the dust opacity law in iHildebrandl (119831 1 
and a dust emissivity index of (3 = 1.5, appropriate for r nassive dense cores ba sed on multi¬ 
wavelength observations at mm and submm wavelengths (iBeuther et al.ll2007al) . we obtained 
gas masses of the objects ranging from 0.5 to 16 Mq, assuming a gas-to-dust mass ratio of 
100. These parameters are also listed in Table 1. The corresponding H2 densities in these 
sources range from 10® to 10^ cm“®. 


3.2. Molecular Line Emission 


Previous arc-second resolution observations detec 


G28.34-P1 region with t he VLA and SMA, respectively (IWang et al.ll2008l: IZhang et al.ll2009 


ed NH3 and emission in the 


Wang et al.ll201ll. 12012 ). Besides CO, the SMA observations in the 230 GHz band did not 
detect line emission at a la rms noise of 100 mJy beam“^ (at a resolution of 1".2 and a 
spectral resolution of 1.3 km s“^). The ALMA observations in the same band reveal line 
emission from 19 molecular species including I80, ^^CS, N2D+, DCN, H2S, SiO 5-4, H2CO, 
CH3OH, and SO, at a level of > 5 ct with a Icr rms noise of 2.5 mJy beam“^. 

Figure 3 presents the line spectra from the SPWO spectral window covering 215.7 
through 217.5 GHz toward Components 9 and 38, the brightest continuum peaks in Core 
2 and Core 5 (see Table 1), respectively. The fluxes are corrected for the ALMA primary 
beam attenuation for direct comparison. Despite the fact that the two cores have a similar 
flux in dust continuum, the molecular line emission from the two are very different. The 
line emission from Core 5 is much fainter than from Core 2. Of particular note is that the 
emission from I80 and SO is nearly non detected in Core 5. A lack of I80 emission is consis¬ 
tent with CO d epletion in co l d and dense environments at the very early stage of evolution 
as suggested in IZhang et al.l (120091) . The fact that we do not detect SO in Core 5 may be 
related to a lack of injection of ice mantles in the gas phase by thermal desorption, i.e., by 
an increase of temperature of dust in this core. The chemistry of SO is initiated by the 
injection of H2S from mantles. H2S is de stroyed by H an d H30’'' to form S and nsS"*", which 
then react with O and OH to form SO (ICharnlevi 119971) . A lack of SO emission in Core 5 
indicates that the protostar(s) embedded are at an even earlier stage of evolution than those 
in Core 2. 





























Figure 4 presents the integrated emission of 18o, ^^CS, N 2 D’^, and H 2 S lines. The line 
emission of these molecules is detected mostly near the dust cores. The I 80 emission is 
strongly detected toward Core 1 through Core 4. There is a faint emission toward Core 
5. On the other hand, emissions of ^^CS and H 2 S are seen toward Core 2 through Core 4. 
The N 2 D’'" emission is detected in the dense ridge along the hlament. However, the emission 
appears to avoid the dust continuum peaks. 


Besides the above mentioned molecules that mainly trace gas around dense cores, H 2 CO, 
CH 3 OH and SO are also present towards Cores 1 through 5 (see Figure 4). Strong emission 
from H 2 CO, CH 3 OH and SO line wings > 3 km s“^ from the cloud systemic velocity is 
detected outside of dust co ntinuum emission. Their spatial coincidence with the outflows in 
CO 3-2 flWang et al.ll 201 ll) and 2-1 (see Section 3.4) indicates their outflow origin. We will 
discuss outflows in the region based on the CO 2-1 and SiO 5-4 data in Section. 3.4. Cores 
2, 3 and 4 are seen in I 80 , ^^CS, N 2 D’^, H 2 S, H 2 CO, CH 3 OH and SO, thus appear to be 
chemically more active than Cores 1 and 5. 


3.3. Core structure and dynamical state 


At a much higher sensitiy ity than the previous high resolution studies of IZhang et ah 


020091) and IWang et al.l fj201l[l . the 1.3mm continuum data from ALMA reveal additional 
fragments. Core 1, which is shown as a single dust peak in the SMA observations, ex¬ 
hibits three additional condensations. These structures are likely the outcome of continuing 
fragmentation upon the formation of cores along the hlament. We investigate the ther¬ 
modynamical properties from the parsec-scale clump to the condensations (~ O.Olpc). We 
measure the size of the clump, cores and condensations using a two dimensional Gaussian 
ht to the dust continuum emission. The htting reports the integrated hux and the size of 
the structure. We then measure the line width at the same spatial sc ales. For the clurnp. we 


use the 1.3mm dust continuum data from the IRAM 3 0m telescope flRathborne et al 


and the NH 3 data from the Effelsberg 100m telescop e Jpillai et ah 


For dense cores, we use the SMA 1.3mm continu um (IZhang et al 


2006 


Wang et al. 


200611 


200811. 


2009 !) . and the NH 3 data 


obtained from the VLA (jWang et al.ll2008l. 120121) . At the scale of condensations, no spectral 
line emission was detected before the ALMA observations. In the ALMA data, a group of 
molecular lines, e.g., SO, H 2 CO and CH 3 OH, is affected by protostellar outflows. Emission 
from I 80 , ^^CS, N 2 D’'' and H 2 S is found toward the dust continuum emission. N 2 D+ is in 
general a reliable tracer of cold and dense gas. However, its emission is not present toward 
all condensations, and does not appear to correlate well with the dust continuum peaks. On 
the other hand, it appears that the I 80 emission coincides with the dust continuum emission 
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well. We measure the line width using the 18o data and found a typical FWHM of 1.7 
km s“^. This value is consistent with the line width of N 2 D’*' in the vicinity of Core 3, the 
strongest emission in the map. 

From the measurements described above, we compute the gas mass from the dust con¬ 
tinuum emission and the virial mass for the clump (size ~ 1 pc), cores (size < 0.1 pc) and 
condensations (size < 0.01 pc). The virial mass is computed following 


Myir = 3k 


G 


Here r is the radius and ay is the line-of-sight velocity dispersion, and G the gravita- 
tional constant, k = is a correction factor dependent on the density prohle p oc 


flMacLaren et al.lll988h . For a constant density, a = 0, k= 5/3, the above equation can be 
rewritten as 

pc km 

Here AVfwhm is the FWHM of the line width along the line of sight. Table 2 lists gas mass, 
line width, radius and virial mass of the cloud features identihed. 

As shown in Table 2 , the line width measured from the NH 3 (1,1) line decreases from 
the clump scale of 0.6pc in diameter to 0.05pc in dense cores, and then increases at scales of 
0.017pc at which the condensations are identih ed. The decrease of line width toward smaller 
spatial scales in this region is hrst reported by IWang et al.l (120081 ) when analyzing the NH 3 
line width obtained from the VLA. They found that the NH 3 line width from the VLA alone 
is systematically smaller than the line width fro m the NH 3 spectra from the combined VLA 
and the Effelsberg 100m telescope. In addition, IWang et al.l (120081) found that the NH 3 line 
width is inversely correlated with the NH 3 huxes. The authors interpret the change of line 
width as indications of turbulence dissipation from the clump to core scales. 

The ALMA observations reveal for the hrst time spectral line information in dust con¬ 
densations identihed with the SMA at 880 pm and ALMA at 1.3mm. The line width derived 
from I 80 in condensations are consistently larger than the NH 3 line width in the cores. Since 
star formation has already taken place in these c ondensations, the line width can be brq ad- 
ened by infall and rotation in the envelope (e.g., IZhang fc Holll997l: IZhang et al.lll998alJbl) . 
The I 80 line width can also be enhanced due to an increase of turbulence thanks to the 
injection of energy from protostellar outhows (see next Section). 

The virial parameter, a, dehned as Myir/Mgas, is smaller than 1 in clump, cores as well 
as condensations for a constant density distribution. Pr evious observation s of G28.34 found a 
density prohle p oc with the power-law index, a=2 (IZhang et aLll2009l: IWang et aLll201l[) 
for detected cores. It is reasonable to assume a similar density prohle in condensations. 
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With a power-law index of 2, the virial masses reported in Table 2 are reduced by a factor 
of 5/3. Applying this correction to the a values in Table 2 leads to a much smaller than 1. 
Our observations indicate that molecular gas is not in a virial equilibrium from the clump 
to condensations. 


The virial analysis presented above does not include magnetic helds, which can in turn 
increase the virial mass. Magnetic helds provide additional pressure in the medium. Their 
contribution to the virial mass for a critical mass-to-hux ratio of gravitational collapse is 
given hj Mb = Here Mb is the contribution to the virial mass from magnetic helds. 

There are no direct measurements of magnetic helds in G28.34 cores. However, recently 
studies of dust polarization in regions harboring high-mass protostellar object s found that 


magnetic helds play an important role in the fragmentation of molecular clumps fjZhang et ah 


helds of order of 1-10 mG ( 

Girart et al 

2009 

MTang et al. 

2009; 

Girart et al. 

2013; 

Oiu et al. 

2013; 

Frau et ah 

2014; 

Oiu et ah 

2014) 

which can increase the virial mass signiflcantly (e.g. 

Frau et ah 

2014 

. Recently, 

Pillai et al. 

(2015 

) reported a magnetic held strength of 0.27 


mG in an IRDG clump Gil.11 using the dust polarization data from JGMT. Assuming the 
same magnetic held hux in the G28.34 PI clump, we hnd the mass for a critical mass-to-hux 
ratio of 1.6 X 10^ Mq. The virial parameter including the ehect of magnetic helds becomes 2 
in the clump. We cannot evaluate virial parameters for the cores detected here since there 
is no direct estimate of magnetic helds in IRDG cores. However, for the sake of argument 
we assume a magnetic held of 0.27 mG, the same as that reported for IRDG Gll.ll, the 
virial parameters derived for cores detected here are still smaller than unity. We therefore 
suggest, a held strength of 2 mG - 5 mG would be required to make virial parameters close 
to 1 for these cores. 


3.4. Molecular outflows 


Molecular outflows are detected using several tracers including GO, SiO, H 2 GO, GH 3 OH, 
and SO. Figure 5 presents GO 2-1 and SiO 5-4 emission in outflows. As shown in Figure 
5, high velocity line emission is detected around Gores 1, 2 , 3, 4 and 5. There appears to 
be more than one outflow originating from Gores 1, 2, 3 and 4, consistent with the presence 
of multiple protostars indicated by multiple continuum peaks. We identify outflow pairs by 
inspecting the channel maps of the GO, SiO, H 2 GO, GH 3 OH, and SO emission. A total of 10 
outflows are identifled in the regi on. Of thern , three outflows (Outflows 3a, 3b and 4c) were 
not seen in the previous study bv lWang et ahl fl201ll) . Table 3 lists the tracers through which 
they are identifled. The GO emission is the most effective tracer for most outflows. However, 
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outflows la and 4b exhibit much stronger emission in SiO, CH3OH and H2CO emission. SiO 
and CH3OH are h eavily depleted with low abundances flO~^^ and 10~^°, respectively) in cold 


and dense regions flMartin-Pintado et al.lll992l: Ijimenez-Serra et al.ll2004ii . Their abundances 
are signihcantly enhanced in protostellar outflows due to shocks triggered by the interaction 
between protostellar wind and the ambient medium, and release the Si and CH3OH molecule 
off the dust grain. The enhanced SiO and CH3OH are seen in both l ow-mass and high-mass 
protostellar outflows (e.g. L1157. IZhang et al.lll995l: lOin et al.ll2007l) . 


Some of the outflows in the region display remarkable collimation. In particular, outflow 
2a associated with Core 2 extends more than 40", or 1 pc in the sky. The full extent is not 
known since it is beyond the FWHM of the primary beam of the ALMA 12m antenna. The 
width of the CO outflow is l".l, yielding a collimation factor, dehned as the ratio between 
the major and the minor axis, of more than 35. This outflow, reaching a velocity ± 50km s“^ 
from the cloud velocity, also has a wide angle component seen in the CO Vlsr velocities 
of 90 - 92 km s“^. The component has an opening angle of 10°. The outflow consists of 
a chain of CO knots appearing symmetrically w ith respect to the outflow origin, which is 
also seen in the CO 3-2 emission with the SMA fjWang et al.ll201lh . The projected spacing 
of the CO knots are 2".5, or 1.2 x 10^ AU. These knots likely correspond to an increase 
in mass ejection related to episodic accretion, as rep orted in both low-mass and high-mass 
protostars (e.g.. lOiu et al.l 120071: lOiu fc Zhang|l2009l) . Assuming a jet velocity of 100 to 500 
km s“^, the time scale associated with the episodic ejection is 5 x 10^ to 10^ years. The 
enhanced accretion/mass ejection can be due to disk instabilities triggered by interactions 
with a companion. 

We derive mass, momentum an d energy in molecu l ar ou tflows identihed in the CO 
following the prescription outlined in IZhang et al.l (120011. l2005l) . We assume that the CO 
emission is optically thin, a nd a CO to ab undance of 10“^. We adopt an excitation 
temperature of 18K following IWang et al.l (120121) . When an outflow is seen in both CO and 
SiO, the outflow parameters are derived using the CO data only. For outflows that are seen 
in SiO only, we use the SiO data to deriye outflow parameters assu ming an optically thi n 
approximation and a relative [SiO]/[H 2 ] abundance of 9.2 x 10“^° (ISanhneza et al.ll2012l) . 
This abundance value of SiO yields consistent outflow masses computed from the SiO and 
the CO emission. The outflow parameters are given in Table 3. 

The mass of the G28.34 outflows, without correcting for the optical depth, ranges from 
0.054 to 2.8 Mq. Besides the outflow la, all other outflows are more massive than 0.15 Mq. 
These values are about one ord er of magnitude grea ter than the typical mass of outflows 
powered by low-mass protostars (iDnnham et al.ll2014jl . and are one order of rn agnitude lower 
than the outflow mass associated with massive protostars (IZhang et al.ll2005h . 
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Protostella r outflows are conn ected to mass ejection during the accretion phase of young 
stellar objects fIShang et ahN2007l) . H igh angular resolution observations provide spatially 
resolved images of massive outflows flSu et ahl |2004 iQiu et al.l 120071: iQiu &: Zhang! 1200911 
associated with individual protostars. The inferred ac cretion rate from massive protostellar 
outflows is typically 10“^ Mq yr“^ f Qiu et ah 2011 1. For G28.34 PI, the accretion rate 
inferred from the outflow rate is around 10“^ Mq yr“^ for outflow 2a. In order to reach a 
star of 10 Mq, it takes 10® years for a constant accretion. This time scale appears to be too 
long. Therefore, it is likely that the accretion rate increases with time as well as with the 
protostellar mass. A variable rate of accretion has been proposed in the theoretical model 
of massive star formation (e.g. iMcKee &: Tad 120021) . and is in agr eement with observations 


that suggest mass infall rates increase with the protostellar mass flZhang||2005l) . 


3.5. Velocity structure along the main filament 


Filaments are dominant structures in the interstellar medium and molecular clouds. 


Recent galactic-wide surveys reinforce this scenario with hndings of spectacular ne 


Molinari et al 


2012 


work 


of hlaments in both nearby and more distant molecular clouds fiChurchwell et al.l 12009 
201o[|. The web-like filaments in massive star forming regions ( Liu et al. 


Galvan-Madrid et al.lboiO . 2013 : Busquet et al. 2013) could be part of the hub-hlament 


structure flMversI l2009h that transports gas and dust to fuel massive star formation at the 
center the web. 

G28.34 PI is embedded in a filamentary infrared dark cloud that spans 6 pc on the 


plane of the sky. We analyze the large scale NH 3 data obtained with the VLA fjWang et al. 


20081 1. Figure 6 presents the peak velocity of the NH 3 (1,1) emission in the region. There 
are complex velocity structures along the major axis of the hlament as well as perpendicular 
to the filaments (near MMIO). Toward clump PI, the line peak velocity varies from 78 to 80 
km s“^. To further examine the gas velocity, we present in Figure 7 the position velocity plot 
of the NH 3 (1,1) emission. The cut of the pv diagram is along the main axis of the hlament 
with a reference position (0", 10") from Gore 3, and a position angle of 44°. We also compute 
the centroid velocity using the main component of the NH 3 (1,1) emission, following the 
formulation of the hrst moment Vc = J F{v)vdv/ f F(v)dv. Here, F(v) is the hux density, 
and V is the line-of-sight velocity of the NH 3 (1,1) emission. We used a hux threshold of 
5 ct noise (10 mJy) when computing the moment data to avoid the contamination from the 
satellite hyperhnes. The right panel of Figure 7 plots the line-of-sight velocity along the 
main axis. We also mark the spatial locations of the dust continuum sources. The hlament 
spans a velocity range from 78 to 82 km s“^, which may arise from large scale gas motions 
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in the filament. 


The position velocity diagram and the NH3 centroid velocity shown in Figure 7 display 
velocity shifts along the main axis of the hlament. There are fou r mm continuum clump s 
along the filament, MM4, MM9, MMIO and MM14, reported in Rathborne et ah ( 2010 1. 
MM4, which is G28.34 PI, and MMIO, MM14 are at velocities around 79 km s“^, while 
MM9 are at velocity around 80 km s“^. There is a large scale velocity gradient of 0.6 
km s“^ over an angular scale of 80" from MMIO, MM14 to MM4. In addition, there is a sea- 
saw velocity pattern with a peak-to-peak amplitude of 0.4 km s“^. Similarly, the filament 
near MM9 shows sea-saw velocity patterns with an amplitude of 1 km s“^. The velocity 
pattern is consistent with gas flows along the hlament toward the clumps and cores (e.g. 
Hacar fc Tafalla 2011 1 


4. Discussion 

4.1. Fragmentation of massive molecular clumps and formation of massive 

cores 


High angular resolution observat i ons of mas s ive IRDC clump s _ in the recent past 


(Wang et al. 

2008; 

Zhang et al.l 

200 

9; 

Swift 

2009; 

Rathborne et al. 

2011 

2011; 

Wang et al. 

2011, 

2012. 

2OI4I 

Tan et al. 

2013; 

Sanhueza et al. 

2OI3I) 


reveal the phys¬ 


ical and chemical state of massive star formation at an early stage of protostar formation. 
Despite typical source distances of kilo-parsecs, interferometric observations achieve suf- 
hcient linear resolution to spatially resolve the global thermal Jeans le ngth in molecular 
clumps (O.lSpc fo r a density of 3 x 10^ cm“^ and a temperature of 15K fjPillai et al.l 12006 


Wang et al.l I2OO8II I. Dense cores revealed in dust continuum emission conta in masses at 


20091: 


Wang et al. 


2011 


Sanhueza et al 


a factor of 10 larger than the global 

Jiermal Jeans mass in the clumd^l ( 

Zhang et al. 

Rathborne et al. 

2OI0I 

Sadavov et al. 

2OI0I 

Zhang & Wang 

2011; 

Csengeri et al. 

2011; 


2 OI 3 I: Wang et al. 2014), but comparable with the turbu¬ 


lent Jeans mass of the clumps; the sound speed in this instance is replaced by the turbulent 
line width. Observations at higher angular resolution suggest that the detected cores them¬ 
selves are sub-fragmented. In IRDC G28.34 and IRDC Gll.ll, for example, dusty cores hrst 
identihed at a resolution of O.OSpc were resolved into several continuum peaks at a spatial 
resolution on the order of 0.01 pc (jWang et al.N201l[ I2014J) . These observations indicate 
that the super Jeans cores continue to fragment and lead to the formation of a group of 


^For cores with masses exceeding the thermal Jeans mass, we refer them as super Jeans cores. 
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stars. However, this need not always be the case. Some super Jeans cores rernain s tar- 
less, but exhibit oscillations in the radial direction(e.g. lAnathpin dika & Franc escoll2013r) . as 
demonstrated in hydrodynamic models for such cores (see also. iKeto &: BurkertI 120141) . 


A similar trend is also observed when ex amining the relation between the mass of frag¬ 
ments and the separation. IWang et al.l (120141) analysed observations of 4 IRDC clumps and 
found that the fragment mass and separation follow the turbulent Jeans fragmentation if the 
sound speed in the Jeans formula is replaced by the turbulent line width. This fragmenta¬ 
tion is significantly different to filaments typically observed in low-mass star-forming clouds. 
Pre stellar cores along these f ilaments are usually separated by a thermal Jeans length (see 


e.g. 


Arzoumanian et al.l 120111) . 


These observations offer direct comparison to theoretical models on massive star and 
cluster formation. Observations hnd that massive stars and clusters form in higher density 
and more turbulent re gions of molecular clouds as compared to their low-mass counterparts. 
Bonnell fc Batel (120021 ) propose in the competitive accretion model that clouds fragment ini¬ 
tially into cores of thermal Jeans mass. These cores subsequently form low-mass protostars 
that accrete the distributed gas from a reservoir of material in the molecular clump. Pro¬ 
totars located near the center of the gravitational potential accrete at a higher accretion 
rate because of a stronger gravitational pull, thus, experience a f aster mass growth. This 
competitive accretion model reproduces the stellar IMF observed (iBonnell et al.l 120041) . 


Alternatively, iMcKee fc TanI (120021) put forward a turbulent core model, in which stars 
form via a monolithic collapse of a massive core. In this model, cores are supported by 
turbulence in a virial equilibrium and have masses much larger than the thermal Jeans 
mass. Protostellar feedback such as heating from the embedded protostars increases the gas 
temperature, and thus, suppress fragmentation. Therefore, cores h arboring massive stars 


usual ly do not sub-fragment significantly and form one or a few stars (iKrnmholz et al.ll2005 
20071) . 


Observations of typical cores found in IRDCs suggest that relatively massive cores could 
possibly sub-fragment to form smaller cores which then spawn protostars. The embedded 
protostars continue to accrete material from their respective surrounding environment. This 
picture is similar to the competitive accretion model in that cores do not acquire all their 
mass before the birth of a protostar, but continue to gain mass during the protostellar 
accretion. However, massive cores contain super Jeans mass, unlike what is assumed in the 
competitive accretion model, but it is close to the turbulent core model. In the meantime, 
massive cores appear to fragment and form a group of stars, in contrast with the picture of 
monolithic collapse. The NH 3 temperature measured in these cores are typically < 20 K, 
indicating that protostellar heating is ineffective in suppressing fragmentation in the core (see 
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also iT^nrigmore et al .1120111) . While magnetic fields can potentially play an important role in 
snppressing fragmentation fjPalau et al.ll2013l. 12014 IZhang et al.ll2014l) . tnrbnle nce provides 


signifi cant snpport in these cores. The NH 3 data from the VLA observations flWang et ah 


20081 ) demonstrate that despite tnrbnlence decay, the line widths measnred at 3" — 4" scales 
appear to be large enongh to snpport the cores in G28.34 PI. Once the protostars are formed, 
ontfiows inject energy in to the immediate snrronnding of protostars, and increase the line 
width as seen in the I 80 line width. This feedback provides additional snppo rt that may 


stop farther fragmentation in the gas at the early stage of a clnster formation fjWang et ah 

2 OI 0 I ). 


One of the key assnmptions in the monolithic collapse model is th at clnster forming 
clnmps are approximately in hydrostatic eqnilibrinm flMcKee fc Tanll2002[) . This is based on 
the observational fact that the time scale for star formation is typically several dynamical 
times (the free-fall time scale of the gas). The virial analysis of the G28.34 PI clnmp and cores 
and condensations within the core reveals that these entities are far from a virial eqnilibrinm. 
The virial parameter, defined as M^r/Mgas, is less than 0.47 from the clnmp to cores and 
condensations (See Table 2). This finding snggests that molecnlar gas in the clnmp is not 
in a virial balanc e dnring star formation. Similar findings of snb-viral paramet ers are also 
reported recently f Pillai et ah 2011 : Li et ah 2013; Tan et ah 2013; Lu et ah 2014). However, 
a key component that is not acconnted for in the virial analysis is magnetic fields. A moderate 
field of 0.5 mG or larger can bring the virial parameter to closer to 1. Althongh magnetic 
fields have not been measnred directly in IRDG cores as present day interferometers lack the 
necessary sensitivity, snch a field streng th is reported in more evolved regions snch as hot 
molecnlar cores (e.g. IZhang et al.l 120141) . As ALMA begins to offer continnnm polarization 
capabilities now and line polarization in the fntnre, it is expected that direct constraints 
on magnetic field strengths in regions at the early stages of clnster formation will become 
available in the next several years. 


4.2. Formation of massive stars through low- to intermediate-mass stages 

G28.34 PI has a Inminosity of 10^ Lq, a gas temperatnre of < 20 K at a spatial scale 
of O.lpc, and a relatively small line width of < 1.7 km s“^ in FWHM. In the same complex, 
G28.34 P2, a molecnlar clnmp of 880 Mq with a Inminosity of 10^ Lrr^ and an N H 3 temperatnre 
of 45 K, has embedded protostellar object(s) aronnd 8 M^ flZhang et al.ll20091) . The relatively 
moderate Inminosity, low gas temperatnre and smaller line widths indicate that G28.34 PI 
is at an earlier evolntionary stage than P2 which already has embedded massive protostars. 
The large reservoir of molecnlar gas in PI demonstrates its potential to form stars in fntnre. 
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Based on the empirical mass-size relation (e.g., iKanffmann fc Pillaill2010l) . G28.34 PI will 
likely bear massive stars when accretion is complete. 

While the comparison of physical properties between G28.34 PI and P2 places PI at an 
evolntionary stage prior to the presence of a massive protostar, the spectral line data from 
A LMA farther constr ains the evolutionary stage of PI. The SMA observations presented 


m 


Zhang et ahl (120091 ) did not detect molecular line emission besides the GO 2-1 line in the 
PI clump at a sensitivity of 100 mJy beam“^ (or 0.4 K), in contrast to th e line emission 


from complex organic molecules in G28 P2 (see also IVa.svnnina et ahl 120141) . The ALMA 
observations detect a plethora of spectral lines from 19 molecules at a la flux sensitivity of 
2.5 mJy beam“^ (0.07K). The line emission includes molecules such as GH 3 OH that has a 
low abun dance of ~ 10 ~^° in de nse and cold environment thanks to its depletion to dust 


grains flJimenez-Serra et al.ll2005l) . Protostellar heating and/or shocks from outflows release 
GH3OH to the gas phase and enhance its chemical abundance. Likewise, SO and OGS 
abundances are enhance d as a consequence of the r elease of ice mantles into the gas phase 
by protostellar heating fjJimenez-Serra et ahl 120121) . Therefore, the detection of complex 
organic and sulfur-bearing molecules serves as a valuable indicator of protostellar heating. 

A detailed chemical modeling of the spectral lines detected in G28.34 PI is beyond the 
scope of the present work and best left for a future article. In order to constrain its chemical 
evolution, we resort to observations of an intermediate-mass protostar (L ~ 10^ Lq) in the 
DR21 hlament obtained with the SMA (Zhang et ah, in preparation). The SMA data consist 
of compact and subcompact conhgurations with a synthesized beam of 3". We convolve the 
SMA data to the same linear resolution as that of the G28.34 spectra in Figure 3, and scale 
the fluxes according to (l/H)^ to account for the dist ance difference betw een the two sources. 
Here, we adopt a distance of D=1.5 kpc for DR21 flGirart et ahl 1201311 . Figure 3 presents 
the template spectra of the intermediate-mass protostar in the DR21 filament after the flux 
scaling. The template spectra reveal line emission from molecules I 80 , GH 3 OH, H 2 GO, HG 3 N 
and SO, indicative of protostellar heating in the core. The line fluxes in the template spectra 
match the observed fluxes in Gore 2 well, but are stronger than the line emission in Gore 
5. This comparison indicates that Gores 2 , 3 and 4 harbors intermediate-mass protostars, 
while Gores 1 and 5 are less active chemically, and may have low-mass protostars embedded 
or are at an earlier stage of evolution, although this may be inconsistent with the dynamical 
age of 10^ yrs of the outflow in Gore 5. 
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4.3. Formation of low-mass stars in a cluster 


The clump G28.34 PI, with a mass of 10^ Mq and an average density of 7 x lO'^ cm“^ 
can potentially form a cluster of stars with a total mass of 100 to 300 Mq, assuming a 10% 
to 30% star formation efficiency. If the stars follow a Salpeter IMF, one expects 75 to 223 
stars in the range of 0.5 to 20 Mq. Among these stars, about 10% of them could possibly 
have stellar masses greater than 10 Mq, and 90% of stars have stellar masses less than 10 
Mq, and 60% of stars have masses between 0.5 and 1 Mq. 


The ALMA observations of G28.34 PI reveal strong dust continuum emission along the 
ridge of the hlament. The spectral line data and the presence of complex organic molecules 
indicate that these cores have embedded intermediate-mass protostars. In addition, some 
of the lower-mass fragments (see Table 1) may harbor low-mass stars. We use the dust 
continuum data presented in Figure 2 and Table 1 to construct the probability function of 
the mass distribution. Figure 8 presents the cumulative mass function (A^(> M) ~ M, 
here N is the number of sources in the mass range > M) for the 38 objects identihed. The 
dashed line denotes the slope of the Salpeter IMF. As shown in Figure 8 , if the cumulative 
mass function follows the shape of the stellar IMF, then there is a large dehcit of lower-mass 
cores in G28.34 PI. In the mass range of 1-2 Mq the number of cores is more than 5 times 
lower than the expected Salpeter slope. 


_ I f low-mass stars in a cluster arise froin thermal Jeans mass fragmentation flLarson 

2005 : Bonnell &: Bate 2002; Palau et ah 2013 ). G28.34 PI has a global thermal Jeans mass 


of 2 Mq, and Jeans length of 0.1 pc. The ALMA continuum image at 1.3mm reaches a la 
rms noise of 0.075 mJy beam~ ^. For an average temperature of 15 K derived from NH 3 
observations flWang et al.ll2012l ) , this flux corresponds to a mass of 0.065 Mq, assuming a 
dust opacity law of iHildebrandl (119831) . a dust emissivity index of 1.5, and a dust to H 2 ratio 
of 1 : 100. The Icr rms noise in the ALMA 1.3mm continuum image is a factor of 30 times 
lower than the global thermal Jeans mass in the molecular clump. The linear resolution 
afforded by ALMA observations is 0.02 pc, much smaller than the Jeans length. Therefore, 
a lack of detection of cores 1-2 Mq is signiflcant. 


In order to examine the limitation of dynamic range in the ALMA image, we perform 
simulated observations using the clean model of the 1.3mm continuum image shown in Figure 
2, derived from clean in GASA. To test the ability of recovering low-mass class 0 protostars 
at the distance of G28.34, we also include a low-mass protocluster NGG 1 333 in the model . 
For NGG 1333, we use the 870 pm continuum data from the JGMT archive ( Kirk et ah 2006 ) 
to derive the sky model at the 1.3mm dust cont inuum. While the 1.1mm continuum data 
for NGG 1333 are also available in the literature (lEnoch et al.ll2006[l . and the data are closer 
to the frequency of the G28.34 continuum from ALMA, the 870 pm data are at a resolution 
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of 19".9, which offer better spatially resolved structure for the sky model. Figure 9 presents 
the 870 /im continuum emi s sion o f NGC 1333. Marked in the figure are class 0 protostars 
identihed bv ISadavov et al.l (1201411 . The masses of these cores are 0.5 to 3 Mq. 


To construct the sky model from the JCMT data, we first scale the 870 /rm fluxes to 
1.3mm by a factor of 3 derived from the comparison of CSO and JCMT data. This scaling 
factor corresponds to a dust emissivity index of approximately 1.5. In addition, we scale 
the fluxes of NGC 1333 by {1/DY to account for the distance difference between the two 
sources. From a source distance of 235 pc for NGC 1333 to that of G28.34, the fluxes are 
reduced by a factor of 417, i.e., (4800/235)^. The flux-scaled image is then deconvolved with 
a Gaussian beam with a FWHM of 19".9 to derive a sky model. The linear scale of 19".9 
at the distance of NGC 1333 is approximately the linear resolution achieved by the ALMA 
observations at the distance of G28.34. Therefore, the JCMT data provide adequate source 
structures that match the ALMA observations. 


We perform simulated observations in casa using the sky model. In the simulation, we 
adopt an array configuration with a total of 29 antennas, a precipitable water vapor PWV = 
1.5mm, and an on-source observing time of 34 min, the same parameters during the G28.34 
observations. We also change the coordinates of the NGC1333 sky model to the pointing 
center of G28.34 so that the UV coverage is identical to the G28.34 data. 


Figure 9 presents the simulated image using the clean components of G28.34 only, and 
the one using the clean components plus the sky model of a low-mass protocluster derived 
from NGC 1333. We reverse the East direction in the RA axis for NGC 1333 to avoid low- 
mass cores overlapping with the dust emission along the ridge in G28.34. As seen in Figure 
9b, the simulated image reproduce the G28.34 image in Figure 2 well. In Figure 9c, the 
low-mass protostellar cores added to the model are also recovered robustly. 


Studies of nearby young clusters flKirk &: Mversll201l[l reveal that massive stars are found 
toward the center of the cluster associated with an enhanced population of low-mass stars 
in their neighborhood. ALMA observations of G28.34 PI reveal extended dust continuum 
along the ridge surrounding the compact sources detected with the SMA at 870 /im. It 
is possible that most low-mass protostellar cores form along the dusty ridge. Indeed, the 
ALMA image at 1.3mm reveals spatially extended emission that is not seen in the SMA 
image at the same wavelength. The difference in the integrated flux between the ALMA 
and SMA images is 80 mJy, which is sufficient to form 100 protostellar cores of 0.8 mJy, if 
assuming 100% fragmentation. 


While the spatially extended emission seen in the ALMA data can give rise to a popu¬ 
lation of low-mass protostars in the vicinity of massive stars, a lack of low-mass protostellar 
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cores outside of dense ridge in the PI region is intriguing. The absence of a distributed 
population of low-mass protostellar cores in the clump suggests several possibilities about 
low-mass star formation in a cluster environment: (1) Low-mass stars only form in the im¬ 
mediate neighborhood of massive stars; (2) The distributed low-mass population does not 
form in situ, but instead forms outside of the clump, and follows the global collapse and 
moves to the center of the cluster; and (3) The distributed low-mass protostellar cores have 
not formed yet in G28.34. 


The first possibility can be discounted since it contradicts with th e fact that mos t 
stellar clusters consist of a population of distributed lower mass stars fe.g.. lQiu et ahll2007l) . 
The second possib ility has been reported in numerical simulations of cluster formation (e.g. 
Smith et ahll2013l) . In simulations of global collapse of molecular clumps, gas fragments and 
forms stars along the streams/hlaments. These stars follow the gas flow and fall toward 
the center of the cluster. A similar o bservational picture was also proposed recently flMvers 


20091: iLiu et al.ll201ll: iLu et al.l 120141) . as massive ptotostars are often found near the center 


of radial gas filaments. In G28.34, the parsec-scale velocity field derived from the NH 3 (1,1) 
line appears to be consistent with gas flowing toward the center of the cluster. The velocity 
of the flow, without correcting for the projection effect is < 0.5 km s“^. If the inflowing gas 
and protostars move at 0.5 km s“^, it takes about 10® yrs for the low-mass protostars to 
travel 0.5pc to the central area of the cluster. While we cannot rule out that this process 
may contribute to the formation of distributed low-mass stars in clusters, its effectiveness in 
the G28.34 clump is questionable. In particular, within a 2pc area of PI, the higher density 
gas is mostly distributed within the clump. It becomes counter intuitive that the lower 
density gas outside of the PI clump would form stars while the higher density gas within 
the clump does not. Therefore, we propose that the third possibility is the likely scenario: 
The outer regions of the PI clump have not fragmented to form low-mass protostellar cores. 
The limit reached in the ALMA observations is 0.2 Mq at a 3 ct level. A lack of detection of 
low-mass cores suggests that in clustered star formation, low-mass cores and stars form at a 
later stage, after the formation of massive stars. 


There h ave been disc ussions regarding when low-mass stars form relative to massive stars 
in a cluster flMversI 1201 ill . Since massive protostars evolve in a relatively short dynamical 
time scale of 10 ® yrs, their strong feedback through radiation and ionization may halt the 
growth of low-mass stars. With this consideration, low-mass stars should form prior to 
massive stars. Our observations of IRDG G28.34 PI and other IRDG clumps indicate that 
massive stars undergo evolutionary stages of low- to intermediate-mass protostars, which 
extends their formation time scale. In addition, it appears that massive stars form first 
in a cluster. Regimes that reach high enough densities become gravitational unstable and 
proceed to star formation. The high density regimes are normally reached hrst toward the 
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center of the clump. That explains why massive stars form at the center of a cluster. Thanks 
to a longer accretion time, protostars who have lower mass initially gain mass over time and 
become massive stars. Formation of low-mass stars in the clump follows as gas continues to 
condense while turbulence dissipates in the clump to reach a super critical state. Protostars 
formed early on in the clump benefit from a longer accretion, thus become massive stars. 
The accretion rates inferred from molecular outflows are typically 10“^ Mq yr“^. Should the 
rate be a constant, it takes more than 10® years to form a 10 Mq star. Th erefore, it is likely 
that mass accretion increases over time or accretion could be episodic(e.g. IStamatellos et ah 


2OI2I). 


5. Conclusion 

We present ALMA observations of a massive IRDC clump G28.34 PI. The clump 
has a mass of 10 ^ Mq at a spatial scale of 0 . 6 pc in diameter, embedded in a long filament 
stretching over 6 pc in the sky. We analyse the ALMA data in conjunction with the (sub)mm 
continuum data from the SMA and the NH 3 data from the VLA to assess physical properties 
from the hlament to the parsec-scale clump, to the O.lpc dense cores and to the O.Olpc scale 
condensations. The main findings are: 

(1) The IRDC hlament exhibits a velocity pattern consistent with mass accretion along 
the hlament toward clumps and dense cores. 

(2) The 1.3mm continuum data from ALMA reveal 5 cores consistent with previous SMA 
observations. These cores are at least a factor of 10 more massive than the thermal Jeans 
mass, implying that turbulence and perhaps magnetic helds are important in supporting 
massive cores during the fragmentation. 

(3) For the hrst time, the ALMA data reveal spectral line emission from 19 molecules 
including ^^CO, I 80 , CH 3 OH, ^^CS, H 2 CO, N 2 D+ and SO in the dense cores. Comparison 
with spectral line data from nearby protostellar cores with embedded intermediate-mass 
stars indicates that C28.34 PI undergoes active massive star formation currently at an 
intermediate-mass stage. 

(4) The superior hux sensitivity in the ALMA continuum data reveals additional frag¬ 
ments with masses of a few Mq. Despite a la mass sensitivity of 0.065 Mq, there is a lack 
of a wide spread low-mass protostellar population expected from thermal Jeans fragmen¬ 
tation. This hnding indicates that low-mass protostars form after massive protostars in a 
protocluster. 
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(5) Cores detected in this cloud have masses at least an order of magnitude larger than 
the thermal Jeans mass and appear to be mostly supported by a turbulent velocity held. This 
study indicates a successive star formation in a protocluster, with cores harboring massive 
stars form hrst as gas in these cores becomes super critical when gravity overcomes the 
internal support. Star formation spreads as gas in other parts of the clump become super 
critical and collapse to form stars. 


We thank P. C. Myers for enlightening discussions. This paper makes use of the ALMA 
data. ALMA is a partnership of ESO (representing its member states), NFS (USA) and NINS 
(Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with 
the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and 
NAOJ. This research is partly supported by the National Science Foundation of China under 
Grant 11328301. K. W. acknowledges support from the ESO fellowship. X. L. acknowledges 
the support of Smithsonian Predoctoral Fellowship. 1. J.-S. acknowledges funding received 
from the People Programme (Marie Curie Actions) of the European Union’s Seventh Frame¬ 
work Programme (FP7/2007-2013) under REA grant agreement PIIF-GA-2011-301538. This 
study used observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 



Table 1:: Physical Parameters of Dense Gas Structure 


Source 

RA 

DEC 

T 

ca 

Peak 

ca 

^int 

Mass 

size^ 

PA 

Reference‘s 


(h;m:s) 

(d:m:s) 

K 

mJy/bm 

mJy 

Mo 

"x" 

(°) 


1 

18:42:51.21 

-04:03:05.2 

22.2 

1.9 

4.2 

2.6 

1.2 X 0.3 

137 


2 

18:42:51.13 

-04:03:06.06 

18.5 

2.0 

10.5 

7.7 

1.8 X 1.2 

60 


3 

18:42:51.26 

-04:03:06.5 

19.5 

4.0 

5.4 

3.7 

0.5 X 0.4 

76 


4 

18:42:51.19 

-04:03:07.3 

17.9 

9.3 

16 

12.1 

0.7 X 0.6 

87 

SMAl 

5 

18:42:50.96 

-04:03:10.1 

22.8 

1.6 

16.6 

9.9 

3.4 X 1.4 

42 

SMA2d 

6 

18:42:50.84 

-04:03:10.5 

18.6 

1.7 

1.7 

1.2 

— 

— 


7 

18:42:50.92 

-04:03:11.3 

24.3 

0.9 

0.9 

0.5 

— 

— 


8 

18:42:51.20 

-04:03:11.3 

12.3 

0.7 

1.4 

1.6 

1.4 X 0.3 

93 


9 

18:42:50.86 

-04:03:11.5 

20.7 

7.6 

9.8 

6.4 

0.5 X 0.1 

20 

SMA2a 

10 

18:42:50.77 

-04:03:11.6 

16.8 

5.3 

20.1 

16.4 

1.5 X 1.0 

49 

SMA2b 

11 

18:42:50.81 

-04:03:12.5 

18.2 

2.2 

2.2 

1.7 

— 

— 

SMA2c 

12 

18:42:50.93 

-04:03:12.9 

22.7 

0.8 

6.3 

3.8 

2.4 X 1.6 

113 


13 

18:42:50.70 

-04:03:12.9 

16.6 

1.2 

1.2 

1.0 

— 

— 


14 

18:42:50.75 

-04:03:13.3 

16.7 

1.1 

4.5 

3.6 

1.6 X 1.0 

114 


15 

18:42:50.73 

-04:03:13.8 

17.4 

0.8 

1.4 

1.1 

0.9 X 0.3 

123 


16 

18:42:50.63 

-04:03:13.9 

15.7 

0.7 

1.2 

1.0 

0.7 X 0.6 

23 


17 

18:42:50.69 

-04:03:15.4 

16.2 

2.6 

16 

13.4 

1.9 X 1.5 

95 


18 

18:42:50.31 

-04:03:15.5 

19.1 

1.0 

1.0 

0.7 

— 

— 


19 

18:42:50.93 

-04:03:16.0 

17.5 

0.4 

1.3 

1.0 

1.4 X 0.7 

50 


20 

18:42:50.58 

-04:03:16.3 

18.3 

8.9 

13.8 

10.2 

0.7 X 0.5 

71 

SMA3 

21 

18:42:50.61 

-04:03:17.0 

18.3 

2.8 

10.1 

7.5 

1.7 X 0.7 

29 


22 

18:42:50.74 

-04:03:17.2 

15.6 

1.1 

2.9 

2.5 

1.5 X 0.5 

73 


23 

18:42:50.51 

-04:03:17.2 

17.2 

1.2 

19 

15 

3.3 X 2.3 

169 


24 

18:42:50.70 

-04:03:18.2 

14.9 

0.8 

1.5 

1.3 

0.9 X 0.4 

72 


25 

18:42:50.53 

-04:03:18.2 

17.1 

0.9 

0.9 

0.7 

— 

— 


26 

18:42:50.80 

-04:03:18.7 

12.5 

1.3 

4.5 

4.9 

1.4 X 1.0 

125 


27 

18:42:50.20 

-04:03:19.1 

17.8 

0.6 

1.7 

1.3 

2.2 X 0.3 

99 
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Table 1:: Physical Parameters of Dense Gas Structure 


Source 

RA 

DEC 

T 

qa 

^ Peak 

qa 

^int 

Mass 

size^ 

PA 

Reference‘s 


(h;m:s) 

(d:m:s) 

K 

mJy/bm 

mJy 

Mo 

n 

x" 

(°) 


28 

18:42:50.28 

-04:03:20.2 

18.9 

5.0 

5.0 

3.6 


— 

— 

SMA4a 

29 

18:42:50.23 

-04:03:20.3 

18.9 

6.5 

19.4 

14 

1.4 

X 0.7 

87 

SMA4b 

30 

18:42:50.32 

-04:03:20.9 

18.5 

3.8 

21.5 

15.8 

2.1 

X 1.2 

12 

SMA4c 

31 

18:42:50.39 

-04:03:21.0 

18.7 

3.5 

10.6 

7.7 

1.5 

X 0.7 

113 


32 

18:42:50.39 

-04:03:22.5 

17.1 

1.7 

10.7 

8.5 

1.8 

X 1.6 

6 


33 

18:42:50.49 

-04:03:22.9 

16.9 

0.8 

1.4 

1.1 

1.0 

X 0.1 

173 


34 

18:42:50.73 

-04:03:22.9 

11.4 

0.5 

0.8 

0.9 

0.9 

X 0.3 

92 


35 

18:42:50.33 

-04:03:23.1 

15.3 

0.6 

0.6 

0.5 


— 

— 


36 

18:42:50.30 

-04:03:23.8 

14.8 

0.6 

3.6 

3.3 

2.0 

X 1.2 

134 


37 

18:42:49.91 

-04:03:25.2 

11.6 

3.1 

10.9 

12.8 

1.7 

X 0.8 

98 


38 

18:42:49.82 

-04:03:25.5 

9.2 

13.9 

18.2 

27 

0.4 

X 0.4 

37 

SMA5 


“^The fluxes reported here are primary beam corrected. The beam size is 0".85 x 0".64 at a positional 
angle of 89°.82. 

represents sources unresolved. 


°Source names used in IZhang et al.l (|2009l) : IWang et al.l (|201lh 















Table 2: Virial Parameters in the Dense Gas 


Name 

(M0) 

Al/“ 
(km s“^) 

r 

(pc) 

(M0) 

a’^ 

Mb 

(M0) 

^total 

Clump G28-P1 

1000.0 

2.67 

0.30 

444 

0.44 

1637.9 

2.08 

Core 1 

28.0 

1.2 

0.030 

6.93 

0.25 

9.79 

0.60 

Core 2 

21.0 

1.5 

0.021 

9.91 

0.47 

8.199 

0.86 

Core 3 

22.0 

0.94 

0.023 

4.28 

0.19 

9.90 

0.64 

Core 4 

43.00 

1.10 

0.028 

7.07 

0.16 

14.44 

0.50 

Core 5 

20.0 

1.70 

0.01 

6.34 

0.32 

2.03 

0.42 

Condensation 4 

12.1 

1.70 

0.0075 

4.57 

0.38 

1.06 

0.47 

Condensation 9 

6.4 

1.70 

0.0026 

1.58 

0.253 

0.126 

0.27 

Condensation 20 

10.2 

1.70 

0.0069 

4.17 

0.41 

0.880 

0.50 

Condensation 28 

3.6 

1.70 

0.0012 

0.705 

0.20 

0.0251 

0.20 

Condensation 38 

27.0 

1.70 

0.0047 

2.82 

0.10 

0.402 

0.12 


“The line width for cores is measured from the NH 3 (1,1) data observed from the VLA (|Wang et alJl2012ll . 
Line widths in condensations are measured from the C^®0 2-1 data in this paper. 




Ma. 


Oitotai = where Mb is the magnetic virial mass. 














Table 3: Outflows and their Physical Parameters® 


Name 

Tracers 

Mass 

Momentum 

Energy 

L 

fC 

'^dyn 

Mout (10-5) 



(Me) 

(M0km s“^) 

M0(km s“^)^ 

pc 

10 ^ yrs 

M0 yr-^ 

la^ 

SiO 

0.054 

0.28 

1.1 

2.2 

6 .8-6.1 

0.78 

lb 

CO,SiO 

0.65 

6.7 

38 

0.49 

16-27 

4.0 

2 a 

CO,SiO 

3.8 

75 

830 

1.0 

11-14 

32 

2 b 

CO,SiO 

0.15 

2.3 

22 

0.17 

5.6 

2.6 

3a 

CO,SiO 

0.15 

0.73 

2.3 

0.16 

13 

1.1 

3b 

CO,SiO 

2.1 

29 

270 

0.90 

4.5-9.5 

16 

3c 

CO,SiO 

1.2 

23 

253 

0.74 

9.2-32 

12 

4a 

CO,SiO 

0.56 

8.3 

77 

0.51 

7.0-10 

7.1 

4b^ 

SiO 

0.23 

1.4 

5.5 

2.3 

4.1-12 

2.2 

5a 

CO,SiO 

0.30 

1.9 

6.4 

0.40 

41 

0.73 


“Parameters are obtained using CO 2-1 data unless stated otherwise. Data are not corrected for the inclination 
angle of outflows. 

^Assuming an SiO relative abundance of 9.2 x 10“^° (jSanhueza et alll2012ri . 

“If both blue- and red-lobes are present, both dynamical time scales are listed. 
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Fig. la.— Three color Spitzer composite image (red/green/blue = 70/24/8 /rm) showing 
the IRDC G28.34+0.06. The circle close PI marks the FWHM held of view observed with 
ALMA. 
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Fig. lb.— The integrated intensity of the NH 3 (1,1) emission flWang et al.l 120081) in white 
solid contonrs overlaid on the l.Smrn continnnm (iRathborne et al.ll2010f) in color scale. The 
NH 3 image is contoured at 10% of the peak (1 Jy beam“^ x km s“^). The thin dotted lines 
indicates the 50% and 100% of the sensitivity level of the 7 pointing mosaic in NH 3 from the 
VLA. The NH 3 data have a resolution of 5" x 3". The thick dashed circle marks the PI held 
that ALMA observed in this work. The color bar on the right-hand side of the plot denotes 
the l.Smm continuum hux scales in Jy beam“^. 
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Fig. 2.— 1.3 mm Continuum image of IRDC G28.34 PI obtained with ALMA. Data are 
plotted in logarithmic scale. The dashed-dotted circle outlines the FWHM of the primary 
beam of the ALMA 12m antenna. The shaded ellipse at the lower left corner of the panel 
marks the synthesized beam. The star symbols and numbers mark the brightest dust con¬ 
tinuum components listed in Table 1. The wedge at the top of the panel displays the color 
scales corresponding to fluxes (Jy beam“^) in the logarithmic scale. Features toward the 
lower right of the image are at the 10% response of the primary beam, thus are not robust 
due to inadequate clean. 
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Fig. 3.— ALMA spectra from Cores 2 and 5 (Components 9 and 38 in Table 1, respectively) 
in G28.34. The flux scale is corrected for primary beam attenuation. Also shown are spectra 
for an intermediate-mass protostar in the DR 21 filament obtained from the SMA, scaled to 
the distance of G28.34 for comparison. 
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Fig. 4a.— Velocity integrated emission (Moment 0) of molecular lines from C^^O, ^^CS, N 2 D’^ 
and H 2 S. The range of integration covers the entire velocity range of the line emission. The 
contours for C^®0, and N 2 D+ emission are plotted in equal increment of 10 mJy km s“^ 
starting from 10 mJy km s“^. The contours for H 2 S emission are plotted in equal increment 
of 6 mJy km s“^ starting from 6 mJy km s“^. The star symbols denote the brightest dust 
continuum peaks in Table 1. 
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Fig. 4b.— Velocity integrated emission (Moment 0) of molecular lines from SO, H 2 CO, and 
CH 3 OH. The range of integration covers the entire velocity range of the line emission. The 
contours for SO emission are plotting in equal increment of 20 mJy km s“^ starting from 20 
mJy km s“^. The contours for H 2 CO, and CH 3 OH emission are plotting in equal increment 
of 30 mJy km s“^ starting from 30 mJy km s“^. The star symbols denote the brightest dust 
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Fig. 5a.— Molecular outflows detected in CO 2-1 and SiO 5-4 in the G28.34 PI region. 
The blue and red-shifted outflows are plotted in blue and red contours, respectively. The 
data shown in color scales are 1.3mm continuum emission also shown in Figure 2. The black 
dashed lines mark the outflows identified in the region. The CO emission is integrated over 
a velocity range from 54 to 72 km s“^ for the blue-shifted emission, and 96 to 112 km s“^ 
for the red-shifted emission. The SiO emission is integrated over the velocity range from 65 
to 76 km s“^ for the blue-shifted emission,and 82 to 96 km s“^ for the red-shifted emission. 
The contours for the CO emission are plotting in equal increment of 0.5 Jy km s“^ starting 
from 0.5 Jy km s“^. The contours for the SiO emission are plotting in equal increment of 
0.02 Jy km s“^ starting from 0.02 Jy km s“^. The star symbols denote the brightest dust 
continuum peaks in Table 1. 
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Fig. 5b.— 
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Fig. 6 .— Shown in color scales is the centroid velocity (moment 1) map in the IR dark 
region of G28.34 derived from the NH 3 (J,K) = (1,1) emission obtained from the VLA. 
The white contour out lines the 1.3mm continuum emission from the IRAM 30m Telescope 
(iRathborne et al.ll2010l) . The scale bar on the right-hand side of the plot denotes velocity in 
k m s~^. The + sym bol marks the hve continuum peaks (SMAl through SMA5) reported in 


Zhang et al.l (1200911 . 
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Fig. 7.— Left: Position velocity (PV) diagram of the NH 3 (1,1) emission. The cut is along 
the main hlament with an offset close to Core 3. Right: The centroid velocity along the 
PV cut. The arrows mark the locations of dust continuum sources. 
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Fig. 8.— Cumulative mass function of the fragments identified in G28.34 PI. The data are 
plotted in logarithmic scales. The dashed line marks the slope of the Salpeter initial mass 
function. 
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Fig. 9a.— a): Dust continuum image at 870 /im obtaine d from JCMT o f the protocluster 
NGC 1333. The star symbol marks protostars reported in ISadavov et al.l fl2014j) . b): Simu¬ 
lated ALMA observations using the clean components of the 1.3mm continuum for G28.34. 
The simulation adopts the observing conditions and the on-source integration time of the 
ALMA project. The star symbol in panels b) and c) marks the brightest 1.3mm continuum 
peaks shown in Figure 2. c): Simulated ALMA observations made from the clean compo¬ 
nents of G28.34 and the model from the NGG 1333 protocluster. The east-west direction in 
the NGG 1333 model image is reversed. 
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Fig. 9b.— 


b): G28.34 Simulated Image 
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Fig. 9c.— 


c): G28.34+NGC1333 
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